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Abstract 


This  report  documents  numerical  simulations  of  the  thermal  environment  within  the  Calde- 
cott Tunnel  during  a gasoline  spill  fire  in  April  1982.  A numerical  fire  model  developed  by 
the  National  Institute  of  Standards  and  Technology  (NIST)  was  used  to  simulate  the  fire.  The 
model  was  validated  for  this  application  using  several  data  sets  from  large-scale  fire  exper- 
iments, including  fire  experiments  conducted  at  a decommissioned  highway  tunnel  in  West 
Virginia,  and  fire  experiments  performed  as  part  of  the  NIST  Investigation  of  the  World  Trade 
Center  Disaster. 

For  the  Caldecott  Tunnel  fire,  the  peak  calculated  gas  temperatures  within  the  tunnel  were 
approximately  1100  °C  (2000  °F),  and  peak  wall  surface  temperatures  were  approximately 
950  °C  (1750  °F).  The  peak  temperatures  were  influenced  by  the  net  effect  of  spalling  and 
ablation  of  the  concrete  walls  and  ceiling,  and  determining  the  sensitivity  of  the  calculations 
to  this  effect  was  a key  objective  of  the  study. 
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1 Introduction 


In  the  early  morning  of  April  7,  1982,  a gasoline  tanker  truck  travelling  westbound  in  the  No.  3 
bore  of  the  Caldecott  Tunnel  overturned  and  caught  fire  as  a result  of  a multi-vehicle  accident.  The 
tunnel  connects  Oakland  and  Walnut  Creek,  California.  The  National  Transportation  Safety  Board 
(NTSB)  and  the  California  Highway  Patrol  (CHP)  conducted  investigations  of  the  accident.  In 
addition  to  a detailed  account  of  the  events  leading  up  to  the  fire,  their  reports  contain  information 
about  the  tunnel,  the  damage  to  the  vehicles,  walls  and  fixtures,  and  various  other  details  [1,2].  As 
part  of  its  on-going  study  of  the  safety  of  waste  transportation  casks,  the  US  Nuclear  Regulatory 
Commission  (NRC)  requested  that  the  Building  and  Fire  Research  Laboratory  of  the  National 
Institute  of  Standards  and  Technology  (NIST)  perform  numerical  simulations  of  the  fire  to  assess 
the  thermal  environment  within  the  tunnel.  The  results  of  the  simulations  were  used  as  boundary 
conditions  for  a detailed  thermo-structural  analysis  of  a transportation  cask  [3]. 

The  simulations  described  in  this  report  were  not  intended  to  replicate  every  detail  of  the  fire, 
since  the  information  about  the  gasoline  spill,  tunnel  lining  materials,  etc.,  was  not  known  to  a high 
enough  level  of  certainty  to  permit  an  exact  reconstruction  of  the  event.  The  approach  taken  was 
to  use  what  information  was  known  about  the  incident  as  a starting  point  for  the  calculations,  and 
then  to  vary  the  unknown  parameters  to  ascertain  the  range  of  possible  outcomes. 

A similar  study  to  that  presented  here  was  conducted  by  NIST  to  estimate  the  thermal  envi- 
ronment of  the  Howard  Street  Tunnel  in  Baltimore,  Maryland,  following  the  derailment  in  July 
2001  of  a freight  train  and  the  burning  of  spilled  tripropylene  and  the  contents  of  surrounding  rail 
cars  [4].  The  similarities  and  differences  between  the  Caldecott  and  Howard  Street  Tunnel  fires  are 
discussed  in  this  report. 


1 


* 'iiM 


V ■^.  ''.■.;>r'--;  ^ ■:  ' -r  ""  ' 


- -'i^^ 


.,  g-i 

aoifoiilm?#!  I 

^ 1^;^  Soim/sb  t Mt>iJffefi 
r^^«SD^b  OilTtu>^ 

'so  ^iiif m^mttonW^  fmjr>?^  ^ 


mh  :Miini,^-m  '^rM  ' joamrp  i»trb  o|  ybj^M  A 

viul  n:  ia:^fit)^fmb  ^if/  m I'lsfriii/T  femwoli  lo  Jnyfwo*! 

imi  gnti'riiif^D'i^  iu  q'c:,yuj:j  ^ifij:m  ^^t';!^/^(0■Ti^;rb.iii>(liC^^h^^  $0im4^i  tftb  bnank^  fl'hr  fijiM 

510  ';►#  bfjji'/r  ra;^='jjbh't»-v^ii  i ‘'jiui  'f^»oi*w,'lo^  «?a3^  ^^llb  b^J  4b|)i;4o:» 

■ ' ’ ^,  .l-st«|0ni  ^idf  ni  b^MmMb 


2 Technical  Approach 

In  cooperation  with  the  fire  protection  engineering  community,  a numerical  fire  model,  Fire  Dy- 
namics Simulator  (FDS),  is  being  developed  at  NIST  to  study  fire  behavior  and  to  evaluate  the 
performance  of  fire  protection  systems  in  buildings.  Version  4 of  FDS  was  publicly  released  in  the 
summer  of  2004.  The  following  is  a brief  description  of  the  major  components  of  the  model.  De- 
tailed information  regarding  the  assumptions  and  governing  equations  associated  with  the  model 
is  provided  in  the  Technical  Reference  Guide  [5]. 

Hydrodynamic  Model  FDS  solves  numerically  a form  of  the  Navier-Stokes  equations  appropri- 
ate for  low-speed,  thermally-driven  flow  with  an  emphasis  on  smoke  and  heat  transport  from 
fires.  The  core  algorithm  is  an  explicit  predictor-corrector  scheme,  second  order  accurate 
in  space  and  time.  Turbulence  is  treated  by  means  of  the  Smagorinsky  form  of  Large  Eddy 
Simulation  (LES).  It  is  possible  to  perform  a Direct  Numerical  Simulation  (DNS)  if  the 
underlying  numerical  grid  is  fine  enough.  LES  is  the  default  mode  of  operation. 

Combustion  Model  Eor  most  applications,  FDS  uses  a mixture  fraction  combustion  model.  The 
mixture  fraction  is  a conserved  scalar  quantity  that  is  defined  as  the  fraction  of  gas  at  a given 
point  in  the  flow  field  that  originated  as  fuel.  The  model  assumes  that  combustion  is  mixing- 
controlled,  and  that  the  reaction  of  fuel  and  oxygen  is  infinitely  fast.  The  mass  fractions  of 
all  of  the  major  reactants  and  products  can  be  derived  from  the  mixture  fraction  by  means 
of  “state  relations,”  empirical  expressions  arrived  at  by  a combination  of  simplified  analysis 
and  measurement. 

Radiation  Transport  Radiative  heat  transfer  is  included  in  the  model  via  the  solution  of  the  radi- 
ation transport  equation  for  a non-scattering  gray  gas.  In  a limited  number  of  cases,  a wide 
band  model  can  be  used  in  place  of  the  gray  gas  model.  The  radiation  equation  is  solved 
using  a technique  similar  to  a finite  volume  method  for  convective  transport,  thus  the  name 
given  to  it  is  the  Finite  Volume  Method  (FVM).  Using  approximately  100  discrete  angles, 
the  finite  volume  solver  requires  about  15  % of  the  total  CPU  time  of  a calculation,  a modest 
cost  given  the  complexity  of  radiation  heat  transfer. 

Geometry  FDS  approximates  the  governing  equations  on  one  or  more  rectilinear  grids.  The  user 
prescribes  rectangular  obstructions  that  are  forced  to  conform  with  the  underlying  grid. 

Boundary  Conditions  All  solid  surfaces  are  assigned  thermal  boundary  conditions,  plus  informa- 
tion about  the  burning  behavior  of  the  material.  Usually,  material  properties  are  stored  in  a 
database  and  invoked  by  name.  Heat  and  mass  transfer  to  and  from  solid  surfaces  is  usually 
handled  with  empirical  correlations. 

The  calculations  performed  for  the  Caldecott  Tunnel  fire  used  a computational  grid  whose 
cells  were  about  0.4  m near  the  fire  source.  A coarse  mesh  was  used  farther  from  the  fire,  where 
it  was  not  as  important  to  capture  the  detailed  mixing  of  fuel  and  oxygen.  The  objective  of  the 
calculations  was  to  estimate  the  temperatures  within  the  tunnel  and  the  heat  flux  to  surrounding 
objects.  The  ability  of  the  FDS  model  to  accurately  predict  the  temperature  and  velocity  of  fire- 
driven  gas  flows  has  been  previously  evaluated  by  conducting  experiments,  both  lab-scale  and 
full-scale,  and  measuring  quantities  of  interest.  Over  the  past  3 years,  FDS  has  been  used  in  the 
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Federal  Building  and  Fire  Safety  Investigation  of  the  World  Trade  Center  Disaster.  As  part  of  this 
effort,  large  scale  validation  experiments  were  performed,  including  fully-engulfing  compartment 
fires  [6].  Gas  temperatures,  heat  fluxes,  and  the  major  gas  species  concentrations  were  predicted 
by  FDS  to  within  10  % of  experimental  measurements. 

As  for  tunnel  fires,  before  FDS  was  used  to  simulate  the  Howard  Street  Tunnel  [4],  it  was 
validated  against  measurements  made  in  a decommissioned  highway  tunnel  in  West  Virginia.  For 
the  two  full-scale  tunnel  fire  tests  that  were  relevant  to  the  Howard  Street  Tunnel  study,  FDS 
simulations  were  within  50  °C  (100  °F)  of  the  peak  measured  values. 

FDS  can  provide  valuable  insight  into  how  a fire  may  have  developed.  The  model,  however, 
is  only  a simulation.  The  model  output  is  dependent  on  a variety  of  input  values  such  as  material 
properties,  timelines,  geometry,  and  ventilation  openings.  Since  exact  information  about  the  fire 
site,  fuel  load,  and  timeline  is  never  known,  estimations  are  incorporated  into  the  model.  Based  on 
experience  gained  during  the  Howard  Street  Tunnel  study,  the  thermal  properties  of  the  concrete 
tunnel  walls  were  identified  as  the  most  important  in  terms  of  predicted  peak  temperatures.  A 
parametric  study  was  conducted  on  the  effect  of  the  thermophysical  degradation  processes  (i.e. 
spalling,  decarbonization  and  ablation)  of  the  concrete  lining  the  tunnel.  Details  are  included  in 
the  next  section. 
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3 Caldecott  Itinnel  Fire  Simulation  Parameters 


The  relevant  details  of  the  incident  that  have  been  used  in  the  present  study  were  included  in  a 
report  by  Tom  Shelton  for  the  California  Highway  Patrol  [2],  and  a report  by  the  National  Trans- 
portation Safety  Board  [1].  These  reports  provided  information  about  the  tunnel  geometry,  fire 
location  and  damage.  Additional  information  and  insights  have  been  obtained  from  analyses  of  the 
incident  by  D.W.  Larson,  R.T.  Reese  and  E.L.  Wilmot  of  Sandia  National  Laboratories  [7],  and 
H.  Ingason  of  SP,  Sweden  [8]. 

The  simulations  of  the  Caldecott  Tunnel  fire  included  only  those  details  that  were  needed  to 
predict  a plausible  time  history  of  the  tunnel  temperatures  for  a time  period  of  three  hours  following 
the  ignition  of  the  fire.  The  simulation  was  not  intended  to  characterize  the  smoke  concentration 
or  toxicity  that  might  have  caused  injury  or  death  to  the  people  in  the  tunnel  at  that  time. 


3.1  Tunnel  Geometry 

It  was  assumed  that  the  Caldecott  Tunnel  is  rectangular  in  cross  section,  10.5  m (34.5  ft)  wide, 

5.2  m (17  ft)  high,  and  1028  m (3371  ft)  long.  For  most  of  this  length,  the  tunnel  is  relatively 
straight  with  a 4.0  % upgrade  from  west  to  east.  The  upgrade  was  included  in  the  simulation 
simply  by  altering  the  direction  of  gravity  appropriately.  The  tunnel  was  assumed  to  have  been 
sealed  except  for  the  east  and  west  portals.  This  was  a conservative  assumption  because  there  was 
a ventilation  system  in  operation,  and  this  might  have  allowed  heat  to  escape,  lowering  near-ceiling 
temperatures. 

The  tunnel  walls  were  made  of  concrete,  with  an  assumed  specific  heat  of  0.88  kJ/(kg  K),  den- 
sity 2100  kg/m^  and  thermal  conductivity  1.4  W/(m  K)  at  ambient  temperature  [9].  No  effort  was 
made  to  model  the  ceramic  tiles  or  grout,  but  modifications  to  the  model  were  made  to  account  for 
the  chemical,  thermal  and  mechanical  degradation  of  the  concrete,  all  of  which  were  accounted  for 
via  an  assumed  “heat  of  ablation”  [10].  The  heat  of  ablation  is  a term  first  used  by  aerospace  engi- 
neers to  describe  the  degradation  of  spacecraft  re-entering  the  earth’s  atmosphere  at  high  speeds. 
The  term  has  also  been  used  by  nuclear  engineers  in  assessing  the  impact  of  core  meltdown  on 
concrete  containment  vessels.  The  heat  of  ablation  is  defined  as  the  amount  of  heat  dissipated 
per  unit  mass  of  material  undergoing  a series  of  erosion  processes  resulting  in  the  removal  of  that 
mass.  It  combines  the  specific  heat  of  the  concrete  plus  the  heats  of  reactions  of  the  various  degra- 
dation processes.  Schneider  [10]  provides  references  to  various  researchers  who  estimate  the  heat 
of  ablation  for  normal  concrete  to  be  in  a range  between  1000  kJ/kg  to  6000  kJ/kg  depending  on 
the  type  of  concrete  and  the  method  of  measurement  or  analysis.  Table  3.1  lists  some  typical  heats 
of  ablation  and  the  temperatures  at  which  concrete  undergoes  the  process. 

In  FDS,  a one-dimensional  heat  conduction  equation  for  the  concrete  temperature,  Ts(x,t),  was 
applied  in  the  direction  v pointing  into  the  solid  (the  point  x = 0 represents  the  surface) 


dr 


-ks^(0,r)  =q'  + q'-m'  AH  A 


(1) 


where  p^,  q and  ks  are  the  density,  specific  heat  and  conductivity  of  the  concrete;  q'^.  is  the  con- 
vective and  q'j.  is  the  (net)  radiative  heat  flux  at  the  surface,  m"  is  the  ablation  (mass  loss)  rate  and 
AHa  is  the  heat  of  ablation.  The  ablation  process  was  applied  at  the  surface,  thus  the  heat  required 
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Table  1:  Typical  values  of  the  heat  of  ablation  of  various  types  of  concrete  [11]. 


Concrete  Type 

Ablation  Temperature  (°C) 

Heat  of  Ablation  (kJ/kg) 

Basaltic 

1080-  1380 

1500-2300 

Limestone  Common  Sand 

1130-  1430 

2300  - 3200 

Limestone 

1230-2280 

2900-5100 

Siliceous 

1130-  1730 

1600-2700 

to  liberate  mass  was  essentially  extracted  from  the  incoming  energy  flux.  The  ablation  rate  was 
assumed  to  be  an  Arrhenius  function  of  the  surface  temperature^: 

m"  = A p5  (2) 

The  parameters  A and  E were  chosen  to  achieve  an  ablation  rate  of  0.05  kg/(m^  s)  at  a surface  tem- 
perature of  1000  °C  (see  Fig  1).  The  pre-exponential  factor  A was  chosen  to  be  0.1  m/s  so  that  the 
mass  loss  rate  would  become  significant  at  about  700  °C,  the  temperature  at  which  decarbonization 
occurs.  The  heat  of  ablation  was  chosen  to  be  2400  kJ/kg. 

The  chosen  values  of  the  parameters  related  to  ablation  were  based  on  preliminary  calculations 
in  which  the  relevant  parameters  were  varied  to  obtain  a mass  loss  rate  consistent  with  observations 
of  “spalling  depth”  made  by  Shelton  [2].  The  term  “spalling”  is  commonly  used  to  describe  the 
physical  deterioration  of  concrete.  However,  in  the  analysis  to  follow,  the  process  of  “ablation” 
included  all  degradation  reactions,  such  as  spalling  (which  begins  near  100  °C,  decarbonization 
(which  begins  near  700  °C),  in  addition  to  actual  physical  ablation  of  material  (which  begins  near 
1000  °C).  Experience  with  Eq.  (2)  as  a simple  model  of  pyrolysis  has  shown  that  regardless  of 
small  changes  in  the  various  parameters,  the  most  important  parameter  is  the  heat  of  ablation. 
Ultimately,  the  model  balances  the  energy  fed  back  to  the  surface  from  the  hot  gases  with  the 
energy  consumed  in  the  ablation  process  and  the  energy  lost  due  to  conduction  into  the  wall.  The 
surface  temperature  naturally  adjusts  itself  until  this  balance  is  achieved. 

3.2  Fire 

As  for  the  fire,  Shelton’s  report  to  the  CHP  estimated  that  8,600  gal  (32  600  L)  of  the  original 
8800  gal  (33  300  L)  of  gasoline  were  unaccounted  for  after  the  fire  and  were  assumed  to  have 
evaporated  and/or  burned.  Experiments  performed  by  Ingason  in  Sweden  [8]  suggested  that  a fire 
from  a leak  in  the  gasoline  truck  would  most  likely  have  melted  the  upper  aluminum  lining  of  the 
tank,  exposing  a large  pool  of  gasoline  to  the  increasing  heat  of  the  fire.  Ingason  estimated  the  pool 
to  have  been  roughly  30  m^  in  area.  Larson  et  al.  at  Sandia  National  Labs  [7]  estimated  the  bum 
time  to  be  40  min.  Assuming  the  gasoline  to  have  a density  of  0.75  kg/L  [9],  the  evaporation  rate 
would  have  been  about  10  kg/s.  Preliminary  simulations  of  the  fire  were  performed,  and  none  of  the 
results  contradicted  the  cited  observations  and  analyses,  except  for  the  fact  that  in  the  simulations 
there  was  only  a brief  period  of  burning  at  the  east  portal,  even  though  Shelton’s  report  suggests 
intermittent  flaming  for  the  duration  of  the  fire. 

Hhis  is  the  same  functional  form  as  the  solid  burning  model  within  FDS.  Here,  however,  the  “activation  energy,” 
E,  has  been  lumped  together  with  the  universal  gas  constant. 
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Figure  1 : Ablation  rate  as  a function  of  concrete  sur- 
face temperature. 

4 Calculation  Results 

Simulations  of  the  Caldecott  Tunnel  fire  were  performed  with  various  mesh  sizes  and  values  of 
the  most  important  physical  parameters.  The  intent  of  these  preliminary  calculations  was  to  assess 
the  sensitivity  of  the  results  to  each  input  parameter.  Of  these  calculations,  two  were  chosen  and 
the  results  were  used  as  boundary  conditions  for  a detailed  thermal  analysis  of  a transportation 
cask  [3].  One  of  the  calculations  can  be  described  as  a “best  estimate,”  that  is,  the  input  parameters 
were  chosen  based  on  the  information  gleaned  from  the  various  post-accident  reports  and  the  most 
likely  thermo-physical  properties  of  the  fuel  and  concrete.  The  second  calculation  was  designed  to 
serve  as  an  upper  bound  in  terms  of  gas  and  surface  temperatures.  Specifically,  it  was  assumed  that 
the  concrete  did  not  undergo  any  spalling  or  ablation.  This  is  not  a realistic  assumption  because 
the  eyewitness  accounts  described  considerable  “spalling”  throughout  large  sections  of  the  tunnel. 
Nevertheless,  it  was  useful  to  provide  a most  plausible  temperature  field,  along  with  an  upper 
bound.  Details  of  the  input  parameters  used  in  the  simulations  are  given  in  Appendix  A. 

Predicted  time  histories  of  gas  and  surface  temperatures  are  provided  in  Figs.  2 and  3.  The 
highest  gas  temperatures  were  predicted  downstream  of  the  fire  near  the  ceiling  (“Upper  Layer 
Gas  Temperature”),  along  with  the  highest  surface  temperatures  (“Ceiling  Surface  Temperature”). 
Temperatures  at  various  other  locations  are  included  for  comparison.  Also  note  the  inclusion  of  the 
ASTM  E 1 19  and  UL  1709  standard  time-temperature  curves  [12,  13].  These  temperature  curves 
are  used  in  fire  test  furnaces  to  evaluate  the  performance  of  building  assemblies.  The  ASTM  E 
1 1 9 curve  has  been  used  for  over  80  years  to  characterize  typical  building  fire  scenarios.  The  UL 
1709,  or  “Rapid-Rise”  curve,  is  a more  recent  addition,  meant  to  characterize  a rapidly  growing 
hydrocarbon  fuel  fire.  Obviously,  it  is  more  representative  of  the  simulated  Caldecott  Tunnel  fire, 
but  it  is  not  intended  to  be  an  upper  bound. 

Some  of  the  major  results  of  the  simulations  were: 

• Both  simulated  fires  (with  and  without  ablation)  induced  a 3.3  m/s  flow  from  west  to  east 
(spatially  and  temporally  averaged).  The  gases  accelerated  downstream  of  the  fire  due  to 
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thermal  expansion.  Predicted  time  histories  of  gas  velocities  at  various  locations  in  the  tunnel 
are  shown  in  Fig.  4. 

• Both  fires  consumed  roughly  70  % of  the  available  oxygen,  with  a heat  release  rate  of  about 
400  MW.  This  is  probably  an  over-estimate  because  the  model  uses  a simple  “mixed  is  burnt” 
combustion  model  in  combination  with  an  empirical  local  extinction  algorithm.  The  actual 
combustion  processes  are  far  more  complicated  and  potentially  much  less  efficient  in  the 
tunnel  environment. 

• The  maximum  ceiling  surface  temperature  was  about  950  °C  for  the  simulation  with  ablation, 
and  1 100  °C  without  ablation. 

• The  maximum  predicted  gas  temperature  near  the  ceiling  was  just  below  1 100  °C  with  abla- 
tion, and  1 150  °C  without.  This  high  temperature  region  was  located  roughly  40  m to  120  m 
east  of  the  overturned  truck,  near  the  ceiling  along  the  tunnel  centerline. 

• The  simulations  predicted  that  flames  exited  the  east  portal,  but  only  for  a short  time.  Shel- 
ton’s report  does  not  say  for  how  long  flames  were  observed  at  the  east  portal.  The  model 
probably  over-predicted  the  combustion  efficiency  of  the  fire,  in  which  case  most  of  the 
fuel  was  consumed  somewhere  in  the  tunnel,  or  never  consumed  at  all.  Another  possibility 
was  that  the  observed  flaming  at  the  east  portal  was  a result  of  unsteady  evaporation  of  the 
gasoline.  It  was  assumed  that  the  gasoline  evaporated  at  a constant  rate  for  40  min  (about 
10  kg/s).  However,  had  there  been  periods  of  greater  evaporation,  this  would  explain  the 
discrepancy  between  the  observations  and  the  simulation. 

Some  discussion  of  maximum  temperature  is  needed.  In  the  numerical  model,  the  gas  tem- 
perature is  literally  the  temperature  of  the  smoke-laden  hot  gases,  and  the  surface  temperature 
is  literally  the  temperature  of  the  solid  wall  at  the  surface.  Shelton’s  report  suggested  that  the 
maximum  gas  temperature  could  not  have  exceeded  the  melting  temperature  of  copper  (1065  °C) 
because  copper  wiring  in  the  upper  wall  light  fixtures  was  not  melted.  Depending  on  their  location 
relative  to  the  walls/ceiling,  the  temperature  of  these  wires  could  have  been  influenced  by  both  the 
gas  temperature  and  the  nearby  surface  temperature.  The  fact  that  the  estimated  wire  temperature 
of  1065  °C  was  between  the  predicted  peak  gas  and  surface  temperatures  is  consistent  with  the 
numerical  predictions.  It  would  be  difficult  to  be  more  precise  in  this  assessment  without  knowing 
more  details  about  the  exact  proximity  of  the  wires  to  the  wall  or  fixture,  and  also  the  extent  of 
conduction  of  heat  along  the  wire  into  the  wall. 
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Figure  2:  Predicted  gas  temperature  time  histories  at  various  locations  inside  the  tunnel  for  simu- 
lations with  and  without  the  inclusion  of  concrete  ablation. 
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Figure  4:  Predicted  gas  velocity  time  histories  at  various  locations  inside  the  tunnel  for  simulations 
with  and  without  the  inclusion  of  concrete  ablation. 
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5 Discussion 


Tunnel  fires  have  been  a subject  of  interest  in  the  fire  research  community  for  decades.  A common 
question  asked  after  each  fire  is:  “How  hot  did  it  get?”  In  this  report,  both  numerical  modeling  and 
physical  evidence  have  been  used  to  answer  that  question  for  the  Caldecott  Tunnel  fire.  But  what 
about  others?  Are  there  similarities? 

Experiments  were  recently  conducted  in  an  abandoned  highway  tunnel  in  Norway  [14].  The 
tunnel  is  1600  m long,  6 m high  and  9 m wide.  The  fires  were  located  560  m from  the  tunnel 
exit,  and  fans  were  used  to  create  a longitudinal  flow  of  3 m/s.  In  one  test,  10900  kg  of  wood 
and  plastic  pallets  fueled  a fire  whose  peak  heat  release  rate  was  measured  to  be  200  MW  and  its 
peak  near-ceiling  temperature  was  1350  °C.  The  tunnel  geometry  and  fire  location  were  similar  to 
the  Caldecott  Tunnel,  and  the  longitudinal  flow  speed  was  similar  to  that  predicted  by  the  model. 
The  obvious  question  to  ask  is:  Why  were  the  predicted  temperatures  in  the  Caldecott  Tunnel 
fire  (1100  °C)  less  than  those  measured  in  Norway  when  the  tunnel  geometry  and  ventilation 
conditions  were  so  similar?  A possible  explanation  is  the  composition  of  the  tunnel  walls  and 
ceiling.  The  Runehamar  Tunnel  in  Norway  was  lined  with  a fire-rated  insulation  material  for  the 
fire  tests,  whereas  the  Caldecott  Tunnel  was  lined  with  normal  concrete  covered  with  ceramic  tile. 
The  severe  spalling  of  the  concrete  observed  by  Shelton  after  the  Caldecott  fire  suggested  that 
the  walls  could  not  have  sustained  the  high  temperatures  measured  in  the  Runehamar  fire  tests. 
Spalling  and  ablation  would  have  continually  exposed  the  fire  to  cooler  surface  temperatures. 

Next,  consider  the  Howard  Street  Tunnel  fire  in  Baltimore,  July  2001  [4].  Using  the  same 
numerical  model,  the  predicted  temperatures  in  the  Howard  Street  Tunnel  were  less  than  those 
predicted  in  the  Caldecott  Tunnel.  Why?  The  issue  in  the  Howard  Street  Tunnel  fire  was  not 
tunnel  lining,  but  rather,  ventilation.  The  Howard  Street  Tunnel  has  a 0.8  % upgrade  in  the  portion 
of  the  tunnel  where  the  fire  occurred,  whereas  the  Caldecott  tunnel  has  a 4 % upgrade  where  the 
fire  occurred.  The  longitudinal  flow  in  the  Caldecott  tunnel  fire  was  one-directional,  as  evidenced 
by  the  lack  of  observed  sooting  on  the  walls  downhill  of  the  fire.  During  the  Howard  Street  fire, 
smoke  was  observed  pouring  out  of  both  portals.  In  the  simulations  of  the  Howard  Street  fire,  the 
make-up  air  from  both  portals  mixed  with  combustion  products  so  that  the  air  reaching  the  fire  had 
a reduced  oxygen  concentration.  It  is  difficult  to  simulate  exactly  the  combustion  processes  in  such 
an  under- ventilated  fire,  but  it  is  clear  that  the  fire  would  not  have  been  able  to  grow  to  a size  that 
was  measured  in  Runehamar  (200  MW),  and  estimated  in  Caldecott  (300  MW  to  400  MW),  both 
of  which  had  one-directional  flow. 

The  two  most  important  factors  in  a tunnel  fire  are  ventilation  and  lining  material.  Fuel  load, 
so  long  as  it  is  plentiful,  is  not  a limiting  factor  on  temperature,  although  it  will  certainly  limit  du- 
ration. The  highest  temperatures  are  achieved  in  fires  with  adequate  ventilation  and  well-insulated 
walls.  The  Runehamar  experiments  had  both  of  these  elements,  plus  truck-loads  of  combustibles. 
The  Caldecott  Tunnel  fire  had  substantial  ventilation  because  of  its  4 % slope,  but  its  walls  were 
made  of  normal  concrete  and  could  not  withstand  temperatures  much  beyond  1000  °C.  Finally, 
the  Howard  Street  Tunnel  fire  had  neither  good  ventilation  nor  “good”  wall  materials,  and  con- 
sequently exhibited  the  lowest  temperatures  of  the  three  fires,  at  least  according  to  the  numerical 
analyses. 
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6 Conclusion 


The  Caldecott  Tunnel  fire  of  April,  1982,  was  modeled  using  the  Fire  Dynamics  Simulator,  a com- 
putational fluid  dynamics  fire  model  developed  by  the  National  Institute  of  Standards  and  Technol- 
ogy. The  objective  of  the  calculations  was  to  quantify  the  peak  gas  and  surface  temperatures  that 
were  likely  reached  during  the  roughly  40  min  bum  and  the  cool  down  period  afterwards. 

The  peak  calculated  gas  temperatures  were  approximately  1 100  °C  and  the  peak  wall  surface 
temperatures  were  about  950  °C.  These  temperatures  were  consistent  with  physical  evidence  col- 
lected at  the  scene,  especially  copper  wiring  in  the  light  fixtures  that  approached,  but  did  not  reach 
the  melting  temperature  of  1065  °C.  It  has  been  assumed  that  objects  near  the  wall  would  have 
exhibited  temperatures  between  the  peak  gas  and  surface  values. 

A sensitivity  study  was  undertaken  to  ensure  that  variations  in  the  physical  parameters  of  the 
model  and  the  accident  scenario  would  not  lead  to  dramatic  changes  in  the  overall  results.  As 
a check  on  the  uncertainty  in  the  parameters  governing  the  ablation  process,  a calculation  was 
performed  in  which  the  concrete  was  not  assumed  to  spall  or  ablate.  The  results  of  this  calculation 
served  as  an  upper  bound  on  the  estimated  peak  gas  and  surface  temperatures. 
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A FDS  Input  Files 

Below  is  a print-out  of  the  electronic  file  containing  the  parameters  used  by  FDS  to  simulate  the 
Caldecott  Tunnel  Fire.  The  version  of  FDS  used  in  the  study  was  4.05  (February  2005).  Details 
can  be  found  in  Ref.  [15]. 


&HEAD  CHID=' Caldecott TITLE=' Caldecott  Tunnel  Fire'  / 


&GRID  IBAR=1200, JBAR=24, KBAR=12  / 

&PDIM  XBAR0=0 . 0,  XBAR=1028 . 0 , YBAR0=-5 . 25 , YBAR=5 . 25, ZBAR0=0 .0, ZBAR=5 . 2 / 
&TRNX  CC=200 . , PC=400 . / 


&TIME  TWFIN=10800.  / 

&MISC  NFRAMES=5400, REACTION=' OCTANE' , GVEC=-0 .392, 0 .0,-9.802, 
SURF_DEFAULT=' CONCRETE' , AUTOMATIC_Z= . FALSE . / 


&SURF  ID=' FIRE' , HRRPUA=14967 . , RAMP_Q=' fireramp' , RGB=1, 0, 0 / 40  min  fire 
&RAMP  ID=' fireramp' , T=  0.0,  F=0 . 0 / 

&RAMP  ID=' fireramp' , T=  60.0,  F=1.0  / 

&RAMP  ID=' fireramp' , T=2400 . 0,  F=1 . 0 / 

&RAMP  ID=' fireramp' , T=2460 . 0,  F=0 . 0 / 


&SURF  ID 

= 

' CONCRETE 

RGB 

= 

0 . 66, 0 . 66 

HEAT_OF_ABLATION 

= 

1500. 

ABLATION_TEMPERATURE 

= 

1000. 

ABLATION_RATE 

= 

0 . 050 

A 

= 

0 . 1 

C_P 

= 

0.88 

DENSITY 

= 

2100. 

KS 

= 

1 . 4 

DELTA 

= 

1.0  / 

I These  4 lines  removed 
I for  no  ablation  case 


&REAC  ID 
FYI 

MW_FUEL 

NU_02 

NU_C02 

NU_H20 

SOOT_YIELD 


' OCTANE' 

'Octane,  C_8  H_18 
114  . 

12.5 

8. 

9. 

0.050  / 


&VENT  XB=510 . 0, 520 . 0, 1 . 0, 4 . 0, 0 . 0, 0 . 0, SURF_ID=' FIRE'  / 

&VENT  PBX=  0 . , SURF_ID=' OPEN'  / 

&VENT  PBX=1028 . , SURF_ID=' OPEN'  / 


&SLCF  PBY=0. 0,QUANTITY=' TEMPERATURE'  ,VECTOR=. TRUE.  / 

&SLCF  PBZ=4 . 8, QUANTITY=' TEMPERATURE' ,VECTOR=. TRUE . / 

&THCP  XYZ=460 . 0, 0 . 0, 5 . 1, QUANTITY^' TEMPERATURE' , LABEL=' u temp  460  m',DTSAM=30.  / 


&THCP  XYZ=480 . 0, 0 . 0, 5 . 1, QUANTITY=' TEMPERATURE' , LABEL=' u temp  480  m'  / 

&THCP  XB=40 . 0, 40 . 0, -5 . 25, 5 .25, 0 . 0, 5 . 2, QUANTITY='MASS  FLOW' , LABEL=' mass  flow'  / 


&BNDF  QUANTITY=' WALL_TEMPERATURE'  / 
&BNDF  QUANTITY=' GAUGE_HEAT_FLUX'  / 
&BNDF  QUANTITY=' BURN_DEPTH'  / 


Indicates  ablation  depth 
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A few  technical  notes: 


• The  parameter  AUTOMAT  I C„Z= . FALSE  . instructs  the  program  not  to  make  any  adjustment 
to  the  basic  mixture  fraction  combustion  model.  There  are  scenarios  in  which  the  fire  is  small 
and/or  the  numerical  mesh  is  coarse,  where  the  stoichiometric  flame  surface  needs  to  be  re- 
scaled. No  such  adjustment  was  needed  for  the  Caldecott  Tunnel  simulations  and  the  feature 
was  turned  off. 

• The  parameter  NFRAMES  controls  the  number  of  times  that  output  data  will  be  dumped 
during  the  calculation. 

• The  4 % grade  of  the  tunnel  was  simulated  by  rotating  the  direction  of  gravity  tan”^  (0.04)  = 

2.3°.  The  new  gravity  vector  was  —9.81  (sin  2.3°, 0,  cos  2.3°). 

• The  ablation  rate  was  governed  by  the  expression 

= A p5  (3) 

The  parameters  ABLATION_RATE  {m"  in  kg/(m^  s),  ABLATION_TEMPERATURE  (T^  in 
°C),  DENSITY  (p5  in  kg/m^),  and  A (m/s)  on  the  CONCRETE  SURF  line  were  used  to  com- 
pute the  remaining  parameter  E.  This  procedure  established  the  functional  dependence  be- 
tween the  concrete  surface  temperature  and  the  ablation  rate.  The  parameter  HEAT_OF_ABLATION 
(kj/kg)  defined  the  energy  loss  per  unit  mass  ablated.  Removing  all  parameters  associated 
with  the  ablation  process  effectively  eliminated  the  ablation  process  altogether  for  the  upper 
bound  simulation. 

• The  HEAT.OF  JVBLATION  was  set  to  1500  kJ/kg  in  the  model.  This  value  did  not  include 

the  energy  required  to  bring  the  concrete  to  the  ABLATION_TEMPERATURE.  The  literature 
values  for  the  heat  of  ablation  given  in  Table  3.1  include  the  sensible  enthalpy  For 

concrete  at  1000  °C,  c^AT  900  kJ/kg. 

• The  RAMP  lines  dictated  that  the  fire  should  grow  to  450  MW  (area  of  FIRE  x HRRPUA) 
in  60  s,  burn  for  40  min,  then  extinguish  due  to  lack  of  fuel.  Actually,  FDS  set  an  evapo- 
ration rate  that  would  ideally  produce  450  MW,  but  the  computed  heat  release  rate  was  less 
(400  MW)  because  of  the  limited  ventilation. 

• The  SOOT-YIELD  was  difficult  to  determine  for  the  fuel  and  the  fire  scenario.  A few  test 
calculations  indicated  that  the  results  were  not  appreciably  sensitive  to  its  value.  Also,  ex- 
periments performed  as  part  of  the  NIST  World  Trade  Center  Investigation  indicated  that  the 
temperature,  rather  than  the  soot  volume  fraction,  of  the  hot  gas  layer  plays  a greater  role  in 
determining  the  heat  flux  to  neighboring  objects  [6]. 

• Several  hundred  virtual  “thermocouples”  were  prescribed  in  the  input  file  to  record  the  time 
histories  of  gas  and  surface  quantities  at  various  locations.  These  have  been  omitted  for 
brevity. 
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